USTC-ICTS-09-01 



Understanding light scalar meson by color-magnetic wavefunction in QCD sum rule 



Kavli Institute for Theoretical Physics China, Key Laboratory of Frontiers in Theoretical Physics, 
Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, P.R.China, 
2 Interdisciplinary Center of Theoretical Studies, USTC, Hefei, Anhui 230026, P.R.China 



In this paper, we study the nonet mesons as tetraquark states with interpolating currents 
induced from the color-magnetic wavefunction. This wavefunction is the eigenfunction of effective 
color-magnetic Hamiltonian with the lowest eigenvalue, meaning that the state depicted by this 
wavefunction is the most stable one and is most probable to be observed in experiments. Our 
approach can be recognized as determining interpolating currents dynamically. We perform an 
OPE calculation up to dimension eight condensates and find that the best QCD sum rule is achived 
when the current induced from the color-magnetic wavefunction is a proper mixture of the tensor 
and pseudoscalar diquark-antidiquark bound states. Compared with previous results, to sigma(600) 
and kappa(800), our results appear better, due to larger pole contribution. The direct instanton 
contribution are also considered, which yields a consistent result with previous OPE results. Finally, 
we also discuss the r\ problem as a possible six-quark state. 

PACS numbers: 12.38.Cy, 12.38.Lg, 12.39. Mk, 12.40.Yx 



In past decades, the question how to validly interpret scalar mesons with their mass below 1 GeV stimulated 
many discussions and controversies In the naive constituent quark model, they are expected to be SU(3)f nonet 
consisting of a quark and an antiquark, with one unit of orbital excitation for positive parity. However, due to the fact 
that the orbital excitation contributes energy about 0.5 GeV, it is difficult to interpret their light mass as well as their 
mass spectrum |2|. Moreover, ao(980) and /o(980) couple to KK channel strongly, which is in contradictory to the 
prediction by naive qq mesons picture. This situation v ery naturally leads to alternative interpretation about these 
mesons, such as tetraquark states d H H, @, & S (M El El El E1E1 El El E3 , which were put forward many 
years ago in [H [2(j. Recently, 't Hooft et. aL p2l{ and authors of |22|] found out new evidence from the instanton 
induced effective Lagrangian, implying that the predominant component of light scalar meson is tetraquark. 

In 1977, using the MIT bag model (23[, Jaffe suggested the existence of a light scalar nonet with masses below 
1 GeV |2C| . This nonet is composed by bound states of diquark and antidiquark. The dominant interaction 
generating the bound state is from one-gluon exchange which induces the following effective Hamiltonian 
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where C > is the strength factor constant, a i and Xi are 2x2 Pauli matrices and 3x3 Gell-Mann color operators for 
the zth quark. This is a simple generalization of the Breit spin-spin interaction to include a similar color-color piece. It 
is also known as "color-magnetic" or "color-spin" interaction of QCD, which was first discussed in the pioneering work 
of De Rujula, Georgi, Glashow (24[. Hereafter, we will call the eigenstates of -He// as color-magnetic eigenstates. The 
eigenf unctions and corresponding eigenvalues of -ff e // for q 2 q 2 system (tetraquark) have been presented in (l9l. [20|. 
In these work, the eigenstate with the largest mass defect is 

|0+, 9) = 0.972 |0+9[1]) + 0.233 |0+9[405]>, (2) 

with 

# e //|0 + ,9) = -43.36C|0 + ,9). (3) 

where + stands for the J p , 9 denotes flavor SU (3) f -nonet, and 9[1] (9[405]) represents the nonet belonging to 
[l]-representation ([405]-representation) of color-spin SU(6)cs- Explicitly, they are 

|0 + 9[1]) = y||(6,3)3; (6,3)3; (1, 1)) + ^|(3, 1)3; (3,1)3; (1,1)), (4) 
|0+9[405]> = ^1(6,3)3; (6,3)3; (1, 1)> - ^1(3, 1)3; (3,1)3; (1,1)). (5) 

In the right hand side of above equations, there is state of | (6, 3)3; (6,3)3; (1, 1)), where (6,3)3 indicates that the 
diquark is in 6-dimension symmetric representation of color SU(3)c with spin S = 1 (so 2S +1 = 3), and in 3- 
dimension 3 representation of flavor SU(3)f. While (6,3)3 means the antidiquark is in the conjugate representation. 
And (1,1) means the bound state of diquark and antidiquark is singlet both in color and spin. In the following, without 
ambiguity, the diquark and antidiquark will be denoted according to their SU(3)c representations. For example, 6 C 
diquark signifies the diquark's wavefunction is (6,3)3. Similarly, |(3, 1)3; (3,1)3; (1,1)) is comprised of spin-0 3 C 
diquark and 3 C antidiquark. 

Basing on Eq. ([3]), Jaffe claimed that the scalar tetraquarks with masses below lGeV exist and the color-spin part 
of their wavefunctions can be described by |0 + 9). Utilizing the latest data, Jaffe's statement could be roughly checked 
for a visual comprehension. For instance, a data fit of charmed baryons determines the constituent quark masses 



ml w m c d w 360MeV m c s » 540MeV, (6) 

where c is the abbreviation of "constituent" . The strength factor constants related to the light quarks are 

C « C g? «20MeV, with qe{u,d}, (7) 
(C qs = 15MeV, C ss = lOMeV). 

Then, if we assume er(600) as one member of + -tetraquark nonet, the mass of <r(600) could be roughly estimated: 



x " ml - CJ2( X * ■ ■ ~ 4 x 360MeV - 43.36 x 20MeV » 573MeV. (8) 



Obviously, Jaffe's claim is reasonable, and the underlying dynamical consideration should be legitimate. Therefore, 
it is interesting to study Jaffe's tetraquark in the framework of QCD sum rule which relates the nonperturbative 
aspects of QCD to the hadronic physics [28|, |2{|. In other words, we will try to obtain a legitimate QCD sum rule for 
tetraquarks in terms of their color-magnetic eigenfunctions. 

QCD Sum Rule (SR) analysis for scalar nonet mesons as tetraquarks has been widely discussed in the literature 
(e.g., see % H, 0, B H, EH El, El, 0, 0, EE EE E3, El)- Since the correlator of tetraquark- type current operator for 
SR has higher energy dimension than that of ordinary baryon-type one, the operator product expansion (OPE) must 
be considered up to higher dimensional operators (condensates) than ordinary baryons. Technically, it has been widely 
accepted that the OPE contributions from condensates of dimensions higher than eight are very small for tetraquarks 
[l^ . To single scalar tetraquark current, it has been shown in pT| that the contributions from the dimension eight 
condensates are unexpectedly large and become dominant in the left hand sum rule. What is worse, their negative 
contributions break down the physical meaning of the left hand sum rule. In order to solve this problem, in [l2l ]. 
the authors demonstrated that the current including equal weight of scalar and pseudoscalar diquark-antidiquarks 
leads to a strong cancelation of the contributions from dimension eight operators in the OPE, and then gives a good 
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sum rule. In [8[, by assuming mixing of single tetraquark currents, the authors performed a SR analysis for low- lying 
+ -mesons as tetraquarks. However, by now, all work on tetraquark SR has not considered a basic question that 
whether the color-spin-flavor structures of the tetraquark-type currents in SR are consistent with the color-magnetic 
hyperfine interaction mechanism on tetraquarks. The aim of this paper is to pursue this question. 

The key point of this paper is that we think the interpolating current used in SR should inherit a color-spin-flavor 
structure from the color-magnetic wavefunction. This means that we treat a current standing for linear combination 
of 3 C -3 C and 6 C -6 C tetraquarks as the SR interpolating current. We emphasize that this combination or mixture of 
3 C -3 C and 6 C -6 C tetraquarks is determined dynamically by Eq. ([3]) without any additional ad hoc assumptions. Due 
to the non-relativistic nature of color-magnetic interaction, it should be aware of that the induced mixture is specific 
to energy scale around lGeV, which is mass scale of mesons we are interested in. In short, our method is based on 
the well established concept that color-magnetic hyperfine interactions play a crucial role in multiquark physics. 

The strategy of the calculation is what follows. At the first step, we will study the properties of the scalar tetraquark 
5?7(3)/-nonet as color-magnetic eigenstate with the largest mass defect in QCD sum rule by OPE expansion. With 
the method presented in Section 2, we construct interpolating currents that can represent the color-magnetic structure 
of tetra quar k. Then utilizing these currents, and following the standard procedure for tetraquark's OPE calculations 
Hi aim EH El, EH, E3, we obtain the contributions from the operators up to dimension eight. Meanwhile, to 
achieve a reliable sum rule, we require that the pole contributions should reach around 50%. Then we obtain a meson 
mass (600 ± 75)MeV. 

In addition, the instanton effects, in other words the topological fluctuations of gluon fields, play an important 
role in the structure of QCD vacuum [3l[ and spectroscopy of multiquark hadrons [H, [H[. So they should be taken 
into account in the SR calculations. Combining the contribution from OPE and instanton, we obtain a mass about 
720 MeV close to previous OPE results. At this stage, a complete sum rule description of + nonet meson has been 
obtained by us, including both the OPE and instanton effects. 

The paper is organized as follows. In Section II, we will deduce the interpolating currents for + tetraquarks from 
their color-magnetic wavef unctions. In Section III, the analytic results of OPE calculation based on previous currents 
will be presented, followed by the numerical results. In Section IV, the single direct instanton contribution will be 
considered. In Section V, we summarize the results briefly and make a speculation on the extension of our method 
to study mesons with 6 quarks (Fermi- Yang meson). In appendix, we will list some necessary formulas of spectral 
functions and correlators. 



II. INTERPOLATING CURRENT FOR JAFFE TETRAQUARK 

Substituting Eqs. (TJJ and ([5]) into @, we obtain the expression of the color-magnetic wavefunction for Jaffe's + 
tetraquark nonet meson as follows 

|0+ 9) = 0.988|(6, 3)3; (6, 3)3; (1, 1)) + 0.1571(3, 1)3; (3, 1)3; (1, 1)). (9) 

The elements for |0 + ,9) are 6 C , 3 C diquarks and 6 C , 3 C anti-diquarks. Generally, the composite operator for a diquark 
with certain structure of color, flavor and spin is 

£ {-i) p >{-i) p '$ T cr#>, (io) 

where {a, b} and {i,j} are color and flavor indices of quarks respectively. Specifically, qtp — u a , q^ = d a , q^ = s a . 
C is the charge conjugation operator, and T is Dirac matrix determined by the spin of the system. (— l) Pc and (— 1) P/ 
reflect the parities of the diquark's color and flavor wavefunctions respectively. As for wavefunctions being symmetric 
in color or flavor, P c = 0, or Pf = 0, and for anti-symmetric ones, P c = 1 or Pf = 1. Notation {a «-> 6}, {i «-> j} 
represent the color and flavor permutations respectively. Since |(6, 3)3; (6,3)3; (1, 1)) signifies that the diquark and 
anti-diquark are symmetric in color and anti-symmetric in flavor, the composite operator of 6 C diquark can be written 
as 

& T CTq^ - & T CTq^ + $ T CTq^ - g« T CP#. (11) 
In the non-relativistic limit of diquark bispinor q T CTq, spin-1 requires that 

T = {a" u , f, 7 " 7 5 } with a» v = 1 -{^ V - 7V)- (12) 
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Then, inserting (fT2f into (jlip. we obtain all possible composite operators for 6 C spin-1 diquark expressed as below, 
Q»>(6) ee _^( ff W T Ca^g« - <#> T Ca^ + 9 « T Ca^> - <tf> T Co*»tf>) 

= ^(gr^g«-^»^ g f), (13) 

= ±(q® T Ci> q ? -<& T C-f$\ (14) 
q£' } (6) = ^(gi 4) r cy'7 8 ff6 W) - 9 w) r OyV?S° + ^ r C7Vg« - g^ ^W 5 ?* ) 

= 0. (15) 
where 1/(2 V2) is a widely adopted normalization. Likewise, the composite operators of 6 C spin-1 antidiquark are 

q£ j) (6) = ±& i ><rv,Ctf) T -qU)v lu ,C$> T ), (16) 

9? } (6) = ■js&Hctf T - g-«7^gf ) T ), (i7) 

Qif (6) = 0. (18) 

Because |(3, 1)3; (3, 1)3; (1, 1)) means that the diquark and antidiquark are anti-symmetric in color, spin and flavor. 
The composite operators for 3 C spin-0 diquarks belonging to representation (3,1)3 of SU(6) CS x SU(3)f are the 
following ones, 

(3) = i=(g« r C 7 5 g« - g« T C 7 5 g«), (19) 
g«»(3) = -L(g« r CgW- g « r Cg«). (20) 
On the other hand, the composite operators of 3 C spin-0 antidiquarks belonging to the conjugate representation are 

0) = ^(#7 5 cgf } T - ei 5 c$ T ), (21) 



- fo) (3) = ±=(&C& T - T ). (22) 



Q P 



For the time being, we can express the composite operators related to |(6, 3)3; (6, 3)3; (1, 1)) as 

T {ij}{lm} _ gfe) (6) Q^ m ) (6) 

= g» T C^ q ^- q ^a^ T + g« T Ca^e^Cq^ T , (23) 
A ii M lm } ^ Qfe) (6) Q^) (6) 

= g« T eW^ ) 7,Cgf T + g » T C-fq^^C^ T , (24) 

where T, A represent "tensor" and "axial vector" respectively. These notations lie with how the diquark and anti- 
diquark operators vary under Lorentz transformation. In terms of Eqs. (|19[) - (|22p . the composite operators corre- 
sponding to |(3, 1)3; (3, 1)3; (1, 1)) are 

s m{im} = q») (3) q(M (3) 

= e abc e ab , c 4: )T Cfq^q^ ) ^Cq i i) T , (25) 
p {i jH l m} s Q fe) (3) Q^) (3) 

= e ofc e o6 ^f r C#?Wc«f, (26) 
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where S, P stand for "scalar" and "pseudo-scalar" respectively. Following Jaffe, {er, /q, <z+, k} are assumed as 
+ -SU(3) f nonet tetraquarks. For er, since its flavor content is {ud}{ud}, by Eqs. and (f2"4")l . the operators 

corresponding to |(6, 3)3; (6, 3)3; (1, 1)) of a are 

Tq = T^ ud}{acl} = u T a Ca^d b d a a^Cul + u^Ca^ d a u a a^CdJ (27) 
A% = A { 6 ud}{ad} =u T a Crd b dal^Cul + ulC^d a u alfl CdJ. (28) 
By Eqs. (j25j) and (|26)l . the operators corresponding to |(3, 1)3; (3, 1)3; (1, 1)) of a are 

SS = St d}{ud} = ea^'ulC^dcU^CdF, , (29) 
PS = Pt d}m = ta^.ulCd^CF,. (30) 
Similarly, for fo, because of its flavor content -^({us}{us} + {ds}{ds}), the results are 

r -4= (^Us\^Us\-\-^ds\^ds\) \ m m rn m 

T/°=T 6 ^ U = -^(^Ca^s^aa^C^ + u^^SaUaa^Csf) 

+ ^={dJCa^s b s a a^CdJ + d T h Co^ s a d a a^Csl), (31) 
v2 
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^C^SbSa^CdJ + dJC^Sada^CsJ), (32) 

c/o c 7f({ tls H a; 0+{d S }{<&}) 1 / T „ 5 _ 5 „_T , T „ 5 j 5^<-Tn /o,\ 

ee = _ £Qbc ( eafc , c , Ub c 7 s c u 6 '7 C*s c , + e QbV d b C* 7 s c d fe /7 Cs c ,), (33) 

Pi° = Sf 2 = —e a Uz a ti c >ulCs c u bl Cs T c ,+e a y cl <%Cs c d y Cs T c ,). (34) 
The results for a+ ({us}{ds}), k ({ud}{ds}) are the following ones, 

T a + _ T ({us}{ds}) = u T Ca ^ s J a(J ^T + u T C(J ^ s J aCJiiv Cs^ (35) 

A a + ^ A ({„ S }{^}) = u T CrSbdal ^T + u T Crs J al ^T^ (36) 

5 a + ^ ^-»^) =eafcc£Qb , c , u T C7 5 ScJb , 7 5 Cs -T j (3?) 

p a + s p ( {us} {^}) = tabc e ab ,, u lc S JyCs T ,. (38) 

T « = T ({ud}{ds}) = u T Cn ^ db - aa ^ cS r + ulCa^daSaa^CdJ, (39) 

A n = = u T C ^ db s alli CdJ + ulCrd a S a ^CdJ, (40) 

sff s = e^^afterfcq, (41) 

p « s p ( { « d}{ ^}) = eabc€ab , c , u T Cdc - SblC jr t (42) 

Subsequently, from above results and basing on Eq. ([5]), we get the desired all possible simplest interpolating 
currents for tetraquark |0 + ,9) as follows 





= <*T X - 


^ps x 


J x 


= "T* - 


^f3P x 


J x 


= aA x - 


^(3S X 


J 2 


= aA x - 


h/3P 3 X 



(43) 

where X can signifies a,K,a+ and /o, with a — 0.988 and (3 — 0.157. We notice that some indispensable contents of 
the best mixed current in [f| disappear here. The reason is that they are forbidden by requiring the wavefunction of 
diquark to be anti-symmetrized (l9l l2fj| . 
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III. QCD SUM RULE ANALYSIS WITHOUT INSTANTON CONTRIBUTION 
A. General formulas for QCD sum rule 

In sum rule analysis, we usually consider two-point correlation functions: 

n{q 2 ) = i f d 4 xe lqx (0\TJ{x)j\0)\0), (44) 



where J is an interpolating current for the tetraquark. We compute Tl(q 2 ) up to certain order in the expansion, which 
is matched with a hadronic parametrization to extract information of hadron properties. At hadron level, we express 
the correlation function in the form of dispersion relation with a spectral function: 

n(<? 2 ) = r p( 2 s) . ds, (45) 

Jo s-q 2 -te 



where 



p(s) = n^Sis- M 2 n )(0\J(x)\n)(n\jH0)\0), 

n 

= 2n f x m s x S(s - M x ) + higher states, (46) 



with the convention 



(0|J(aO|Si) = V2fimi. (47) 
The sum rule analysis is then performed after Borel transforming both sides of Eqs. (|44|) and (I45() . 

U^(M 2 B ) = B Ml n(q 2 ) = J™ e- s / M *p(s)ds. (48) 

Usually, evaluating p(s) by OPE or some other methods, then from Eq. (j4"8|) , one obtains the left hand sum rule 
(LHS). On the other hand, inserting Eq. (|46|) into Eq. (|48]) . one derives the right hand sum rule (RHS). By definition, 

n RHS (Ml) = 2*ftm* x e- m *' Ma » . (49) 
The LHS and RHS are supposed to be equal, so we obtain 

So 

e~ s l M *p{s)d s = 2Trf x m x e- m x/ M B. (50) 



o 

In above expressions, we have chosen a finite threshold So to exclude the contribution from the continuum. Differen- 
tiating Eq. (|5TJ)) with respect to j^-, and dividing it by Eq. (fBTJ)) . finally we obtain the physical mass 

B 

rSo s ,mI ,s d 

Ml = 51 

J °e- s / M Bp(s)ds 

In the following, we study both Eqs. (|4^|) and (|?T|) as functions of Borel mass Mb and threshold Sq. 

B. OPE calculation for + nonet as Jaffe tetraquark 

The cr-correlator can be expressed as follows: 

W{q 2 ) = i ( d 4 xe lq - x (0\TJ a (x)J a H0)\0) 



= a 2 W A ^ + f3 2 U B »P + amts + n%T0. (52) 

where J CT = aA + (3B represents any one of the four possible currents in Eq. (|4"5|) , A represents the composite operator 
related to 6 C -6 C , and B is that associated with 3 C -3 C . is the correlator between A-type content and B-type 
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content. In this section, we will first compute the spectral functions for the correlators through OPE expansion, 
then insert these results into the Eq. (|48[) to obtain the Borel transformed correlators. In the process of calculating 
OPE, we use the following propagators for quarks 0, which contain all the necessary terms for computing tetraquark 
spectral functions. 



(01^(^(0)] |0) 
iS ab . i A" 



27T 2 5 



tX + 



i8 ab m 
48" 



32tt 2 
q -(qq)x^ 



2 S c G^(a" 
i5 ab ™ 2 



1~2 



8ir 2 x 2 



■x with q G {u, d}. 



d ab x 2 
192 



c qcrGq) 



S ab r 



(53) 



iSf(x) = 



(0\T[s a (x)s b (0)}\0) 
iS ab . 



2ir 2 x 4 ' 
i5' 



■ \n -i sab 

^^- 2 ^x + xan-^s) 



ab f 



(ss)x + 



iS ab m 2 
8ir 2 x 2 ' 



S ab x 2 
192 ' 



tGs 



S ab r 



4ir 2 x 2 



(54) 



Actually, OPE computation for tetraquarks is rather long, but it can be performed analytically. A convenient 
formulation for performing this calculation has been presented in @, Q. The MATHMATICA with FEYNCALC 
[40| may be helpful for computation. In the following, we use the notations and formulations in @, We have 
performed the OPE calculation for spectral functions up to dimension eight, which is up to the constant (s°) term 
of p(s). During the calculations, we have assumed the vacuum is saturated for higher dimension operators, such as 
(0\qqqq\0) ~ (0|qg|0) 2 . After finishing the OPE calculation, we obtain the following results for a meson, 



<xOPE 

Pt,t 



1280 

Qm 8 
_I_C__£. 
1 4tt 6 



16tt 6 



.3 , r 21 K 



(qq)m q , 1\{£GG). 



2tt 4 



768 



18(w)» 



ll{g>GG)rr£ 3(gqaGq)r 



64tt 6 



9m q 



15(qq)m 



2tt 6 
30(qq) 2 m 2 



+ 



ll(g 2 GG)m 



64tt 6 



\\(f?GG){qq)m q . 
48tt 4 



(55) 



aOPE 

Ps,s 



2 3 

m q s 



61440tt 6 15367r e 

M! )S + ( M 



m q (qq) 
256tt 6 96tt 4 



(g^GG) 2 
6144tt 6 ' 



,3ro° 
•64tt 6 



12tt 



k 64tt 6 



512tt 6 



16tt 4 



24tt 2 



(g^GG) 



1117, 



1024tt 6 
2 GG)(gg)m, 
384tt 4 



(gqaGq)m q 

32tt 4 
(gqaGq) (qq) , 

12tt 2 



(56) 



ctOPE 

Pt,s 



ctOPE 

Pp,p 



ctOPE 

Pt,p 



<tOPE 
PS,T 



614407T 6 



GG) 2 



3(^GG)m2 S^GG)™ 4 (g>GG)(qq)m q 



1024tt 6 



m 2 s 3 11m 4 
+ ( 



256?r 6 

m q (qq) 



256tt 6 
g 2 GG) 



512vr 



256tt 6 32tt 4 



6144tt 6 



)* 2 



9m^ 
^64^ 



64tt 4 

5(qq)m 



{gqaGq)i 



W) 



->+( 



3m® 



3(gg)r 



32?r 4 12tt 2/ v 64tt 6 

tfGG)(qq)m q (gqaGq)(qq), 



4tt 4 



(g^GG)rr 

5127T 6 



8tt 4 
3{gqaGq)m 



1287T 4 



12tt 2 



ctOPE 

Pp,t 



(g 2 GG) 2 , 3(^00)7 



/ 3(g 2 GG)m 



512tt 6 



1287T 6 



128tt 6 



+ 



16tt 4 



(g 2 GG)(gg)m g , 
32tt 4 



3(g 2 GG)m 2 

1024tt 6 
31m 2 (gg) 2 

24tt 2 



(57) 



(58) 
(59) 



CT OPE 

Pa,a 



7680tt 6 



9 3 
m q s 

1287T 6 

3m 8 



+ ( 



5m 4 
32^6 



(g 2 GG) 
3072tt 6 



9m| 



5( W )m3 9(g 2 GG)m 2 (gqaGq)m q 



16tt 6 



4tt 4 



3tt 2 



)*+( 



3(gg)m5 , 5(g 2 GG)m 4 



8tT 6 



2tt 4 



256tt 6 



7m 2 (qq)' 
3tt 2 



512tt 6 
{f?GG){qq)m q 
64tt 4 



(gqaGq) (qq) 



3tt 2 



), (60) 



aOPE 

Pa,s 

<tOPE 

Pa,p 



ffOPE _ 3(g 2 GG)m 2 (g*GG)(qq)m q 

PS,A ■ s 



1024tt 6 



64tt 4 



<tOPE _ n 

Pp,a — u - 



(61) 
(62) 
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In above equations, (qq) is a dimension d — 3 quark condensate; (g 2 GG) is a dimension d — 4 gluon condensate; 
(gqaGq) is a dimension d = 5 mixed condensate; the strong coupling constant takes its value at energy scale about 1 
GeV, that is the energy scale we are interested in. Long distance bulk properties of physical vacuum are effectively 
parameterized in these vacuum expectation values. At present, according to Eq. (43), we can make use of above 
spectral functions to generate correlator of each kind interpolating current belonging to a. These correlators will be 
the starting point of numerical calculation in the next section. 

In order to prevent the long listing of formulas for spectral functions from obscuring the conceptual content, we 
will put the necessary spectral functions of k, a + and fo into the appendix. 

C. Numerical analysis of QCD sum rule for OPE contribution 

For numerical calculations, we use the following values of condensates [34l. l35l. l36l l37l l38l l39l |42| : 

(qq) = -(0.240 GeV) 3 , 

(ss) = -(0.8 ± 0.1) x (0.240 GeV) 3 , 

(g 2 s GG) = (0.48 ± 0.14) GeV 4 , 

m u = m d = m q = 0.1 x 2.4~ 3 GeV, 

m s (l GeV) = 125 ±20 MeV, (63) 
(g s qaGq) = -M 2 x (qq) , 
M 2 = (0.8 ± 0.2) GeV 2 . 

Figure 1 shows the LHS of four possible interpolating currents of the a meson, as a function of Borel mass squared, 
in the case of infinite threshold. From the definition of Eq. (|48|) . the LHS should be positive quantities. However, 
in practical calculations, the positivity may not be necessarily realized due to the insufficient convergence of OPE 
calculations. In our case, from Figure. [T] we see that current Jf and current J% show better convergence than current 
Jg and current J". 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Borel Mass 2 [GeV 2 ] Borel Mass 2 [GeV 2 ] 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Borel Mass 2 [GeV 2 ] Borel Mass 2 [GeV 2 ] 



Figure 1: LHS of four interpolating currents of a meson, as a functions of Borel mass squared, with so=infinity, in units of 
GeV 10 . 
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To find the current with the best convergence, we have to refer to their Borel transformed correlators in numerical 
expressions, which are: 



j-j-CT(all) 


= 1.9 X 10" 


■ 5 M B ° 


- 1.9 x 10" 


*M B - 


K9.5 x 10" 


Hl% 


+ 3.7 x 10" 


8 K 


- 8.5 x 10" 


8 M 2 


n CT(all) 
iL 2 


= 1.9 x 10" 


' 5 M B ° 


- 2.0 x 10" 


- 8 M* - 


h9.5 x 10" 


Hi% 


- 4.2 x 10" 




- 2.1 x 10" 


8 M 2 


u u(a,ll) 
li 3 


= 3.2 x 10" 


6 M B ° 


- 2.5 x 10" 


°M B - 


I- 1.6 x 10" 


*M% 


- 6.2 x 10" 


6 M B 


- 5.1 x 10" 


6 M%, 


11 4 


= 3.2 x 10" 


- 6 Mi° 


- 2.5 x 10" 


- 9 M| - 


h 1.6 x 10" 


6 M B 


+ 6.2 x 10" 


6 M B 


- 5.1 x 10" 


6 M%. 



(64) 

From these expressions, it is obvious that current shows the best convergence behavior, so we will utilize current 
to compute the physical mass of a. We first choose an infinite threshold to estimate the mass as the traditional 
sum rule has done [29| . In Figure [H we exhibit the behavior of the mass of a meson as the function of Mb for infinite 
and finite sq. In traditional sum rule, if the mass as a function of Mb, has a wide minimum, then the minimum value 
of mass function can be perceived as the real mass of the state. From Figure we observed that M a as a function 



of M B indeed has a minimum with M a ^ m ^ 

tW 



0.59 GeV at M% 



0.079 GeV . At this value of Borel mass, the 



correlation function n^ 11 * 1 = 3 x 10~ 9 GeV iu , so the positivity of LHS is kept. Although M a ^ min j is very close to the 
experimental center value {M a ) ~ 0.6 GeV, the minimum is not wide enough as required. Therefore, to obtain an 
acceptable result, we have to adopt finite thresholds scheme @, B H, EH to repeat the process of computing mass. The 
results for some values of threshold are presented in the right part of Figure 2. We notice that when the mass becomes 
weakly dependent on Mb, the value of mass is around 0.6 GeV. But we also find that as the threshold increases, the 
mass will increase too. This may be due to the fact that a is a broad resonance state. So there must be some criteria 
to help us dictate which value of mass is the most believable one. Combining the points of view adopted bydGJ, 
on judging when an acceptable sum rule is arrived, we postulate the following criteria. 




0.4 0.6 

Borel Mass 2 [GeV : 



Borel Mass 2 [GeV 2 ] 



Figure 2: Mass of a is illustrated as function of Borel Mass squared. The left figure is in the case of infinite threshold, while 
the right one is in cases of finite thresholds. The results corresponding to so =0.5, 0.6, 0.7 GeV 2 are represented by a solid 
line, a dashed line and a dot-dashed line respectively. 

1. The Borel transformed correlation function Ii(M B ) should show a good positivity for almost all values of Borel 
mass. This is usually related the convergence of LHS. 

2. The physical mass should depend weakly on the value of Borel mass in a wide region. In other words, there 
should be a Borel window. 

3. OPE convergence. This is a strong constraint to the lower bound of the M B region. OPE series converge better 
for higher values of M B , so that requiring a good convergence sets a lower limit to M B . To current , we find such 
a lower limit of M B in the following. We first rewrite the spectral function corresponding to as, 

p(° PE) = ^=oC (8 - 2n) s" = S*=oP n > (65) 

where c ( 8 ^ 2 ™) denotes the operators of mass dimension (8 — 2n), p n = cS 8 ~ 2n )s n . From Eqs. (55)-(62), 
we learn that terms p( 3 ' 4 ) are perturbative contributions denoted as p<-P ert \ in other words, they do not con- 
tain condensate. Remaining terms represent contributions from operators of dimension 4, 6 and 8. These 
terms are dominated by condensates including the non-perturbative effect, denoted by p^ 2 \ p^ respec- 

tively. In Fig. [3j we present the relative contribution of p^ 2 \ p^ x \ p^ to the total spectral function 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.i 



BorelMass 2 [GeV ] 



Figure 3: convergence of OPE series of spectral function related to current J2 for so = 0.6GeV . 
pi OPE) . The thick line denotes [J °' 6 (p (pert) + pW)e~ s / M Bds/ f^ 6 p( OP ^e~ s / M B ds], the dashed line signifies 

[J°-\ p (pert) + p (2) + p W) e -s/Ml ds/ JO.G p (OPE) e -s/Ml ds ^ ^ daghed doted ^ represents [ft* (p<*«*) + p^ + 

p {i) + p (°)) e - s / M Bds/ J° 6 p(° PE h- s / M Bds=l}. We see that, for M% > 0.2GeV 2 , the addition of a subsequent term in 
expansion (|65p . brings the curve closer to an asymptotic value (which is normalized to 1). Furthermore, the changes 
in this curve become smaller with increasing dimension. Thus, for so = 0.6GeV 2 , the convergence is satisfied by 
M| > 0.2GcV 2 . For s = 0.5, 0.7, 0.8GeV 2 , convergence limits Af| > 0.2, 0.3, 0.4GeV 2 , respectively. 

4. For a given threshold, the pole contribution should be sufficient large. By choosing suitable Borel mass, this 
can be satisfied. Since the Borel transformation suppresses the contributions from so > Mj|, small value of Af| are 
preferred to suppress the continuum contributions. But Af| cannot be arbitrarily small, or it will spoil previous three 
requirements. To a, we have found such optimal values of Mg for different thresholds. We list the corresponding pole 
contributions in Table I. The pole contribution is defined as 



Pole contribution : 



Jo 



s o -s/Ml 



p{s)ds 



e~ s / M B p(s)ds 



(66) 



Table I: Pole contributions of various threshold. 



so (GeV 2 ) 


0.5 


0.6 


0.7 


0.8 


Mi (GeV 2 ) 


0.2 


0.2 


0.3 


0.4 


Pole (%) 


40 


52 


35 


25 


Ma (GeV) 


0.6 


0.6 


0.7 


0.75 



From this table, we can extract following information that when threshold changes from 0.5 GeV 2 to 0.8 GeV 2 , the 
pole contribution will vary from 40% to 25% correspondingly, but reaches its maximum 52% at M|=0.2 GeV 2 , when 
so = 0.6GeV 2 . That the pole contribution reaches 52% implies that a good sum rule has been obtained. We get 

772 ,j = (600 ± 75)MeV, with Pole contribution(52%), (67) 

where (±75) MeV originates from the error of condensates (see Eq. l6li|) . It is remarkable that the Pole contribution 
is larger than that given in [8J, where the Pole contribution is below 30%. 

Applying the same analysis to meson k, the LHS of four possible interpolating currents of k can be found in Figure 
21 with threshold value sq being infinity. The corresponding numerical expressions are listed below: 



K(all) 
(all) 



nZ Kau > = 1.9 x 10~ 5 M^° - 1.2 x 10" 6 Af| + 6.7 x 10~ 6 M|, - 1.3 x IQ- J M% - 1.2 x 10" 7 M|, 



1.9 x 10" 5 M^° - 1.2 x 10~ 6 M| + 6.7 x lO^A/J, - 1.9 x lir 7 M% - 5.7 x 10" 



IT 



n «(all) = 3 2 x 10 -6 M 10 

K(all) - 3.2 x 10~ 6 M|° 



1.9 x 10" 7 Af^ + 1.7 x IQ- & M% 



n 



5.2 x IQ- & M% - 4.6 x 10~ 6 M 
1.9 x 10~ 7 M| + 1.7 x 10- 6 A/J, + 5.2 x IQ- & M% - 4.6 x 10~ 6 Af 



B ' 

2^ 



(68) 
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Figure 4: LHS of k meson as functions of Borel mass squared with so=infinity in units of GeV 



From Figure [4] and above expressions, we notice that current J-f, which is a proper mixture between tensor and 
pseudoscalar contents, is the best interpolating current. By setting the threshold to be infinity, we obtain an estimation 
for the mass of k. As shown in Figure M K as a function of Mg has a minimum with M re ( m j n ) 



M| = 0.2 GeV 2 . 



At this value of Borel mass, the correlation function H 



1.6 x 10" 7 GeV 



10 



= 0.90 GeV at 
the positivity of 



LHS is also retained. But the minimum is still not wide enough, then the finite threshold analysis should be performed. 
The results are shown in the right part of Figure At the Borel window, the mass of k is close to 0.8 GeV. To 



8 3 '° 

II 

*2 2.5 

M 

■S 2.0 

o 

1 1.5 
= 

* 1.0 
o 

o.o 



0.2 



0.4 0.6 0.8 

Borel Mass 2 [GeV 2 ] 




Borel Mass 2 [GeV 2 



Figure 5: Mass of k is illustrated as function of Borel Mass squared. The left figure is in the case of infinite threshold, while 
the right one is in cases of finite thresholds. The results corresponding to so =0.7, 0.8, 0.9 GeV 2 are represented by a solid 
fine, a dashed line and a dot-dashed line respectively. 

find the best sum rule, following the previous criteria, we find that to k, the convergence limits M% > 0.25GeV 2 for 
,s = 0.8, 0.9GeV 2 and M 2 > 0.225, 0.3GeV 2 for s = 0.7, 1.2GeV 2 , respectively For instance, to s = 0.9GeV 2 , 
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the convergence of OPE series is shown in Fig. [6] 




1.0 ^- ---------- ------------------------- 

0.2 0.4 0.6 0.8 1.0 

BorelMass 2 [GeV ] 



Figure 6: convergence of OPE series of spectral function related to current J2 for so = 0.9GeV 2 . 
The pole contributions for several values of threshold are listed in Table II. When sq — 0.9GeV 2 , M| = 0.25GeV 2 , 

Table II: Pole contributions of various threshold. 



so (GeV 2 ) 


0.7 


0.8 


0.9 


1.2 


M 2 B (GeV 2 ) 


0.225 


0.25 


0.25 


0.5 


Pole (%) 


43 


47 


56 


27 


M K (GeV) 


0.75 


0.8 


0.82 


0.95 



we get a pole contribution 56%. Such a large pole contribution suggests that a good sum rule has been obtained. We 
get the mass of k, 

m K = (820 ± 80)MeV, with Pole contribution(56%). (69) 

This pole contribution is also larger than that given by Q, where the pole contribution approaches 45%. 

Lastly, for a + and fo that are degenerate in OPE calculations, the LHS of four possible interpolating currents are 
shown in Fig. [7J with threshold value so being infinity. Their numerical expressions are the following ones: 



jjO+,/0(all) 


= 1.9 x 10" 


■ 5 Ml° 


- 2.3 x 10" 


6 M|- 


1-4.0 x 10" 


e M e B 


- 5.8 x 10" 




+ 7.2 x 10" 


7 Ml, 




n a+,/0(all) 


= 1.9 x 10" 


- 5 Ml° 


- 2.3 x 10" 


6 M|- 


(-4.0 x 10" 


*M% 


- 1.1 x 10" 




+ 7.7x 10" 


7 Ml, 




n-o+./OCall) 


= 3.2 x 10" 


- 6 Ml Q 


- 3.7 x 10" 


7 Ml- 


h 1.8 x 10" 


6 M 6 B 


+ 4.2 x 10" 


*M% 


- 4.1 x 10" 


6 Ml, 




jjO+,/0(aU) 


= 3.2 x 10" 


■ e Mh 


-3.7 x 10" 


7 M%- 


h 1.8 x 10" 


6 M 6 B 


+ 4.2 x 10" 


- 6 M| 


- 4.0 x 10" 


6 Ml. 


(70) 



From Fig. [7| and above expressions, current J% + seems to be the best one. But when applying the traditional sum 
rule method to estimate mass, it turns out that there is no minimum as shown in Fig. [5] Furthermore, if we choose 
certain threshold and Borel mass to reproduce the experimental center value of the masses of a+ and fo, the pole 
contribution can only be around 10%. This indicates that in contrast to the success of SR analysis of a and k, the SR 
fails to analyze a+ and fo, in terms of the interpolating currents deduced from their wavefunctions as tetraquarks. 
The reason is as follows. Jaffe's wavefunctions are the eigenfunctions of H e ff in Eq. ([TJ) . However, He ft is only 
an approximate description of color-magnetic interactions Hcm = — Yli j • \j)(~<Ji ■ ~a j) @, HI, [26|, [27[. If 

the flavor 5'?7(3)/-symmetry is exact, the interaction strengthes Cij are flavor-(ij') independent, i.e., CV, = C, then 
Hcm = H e ff. But for real QCD, the constituent mass m c u « m c d , while m c s > m c = (m° + m c d )/2. So SU(3)f must 
be broken within order 0{{m c s — m c )/m c s ) ~ O(0.3). Therefore, both H e ff and Jaffe's wavefunction |0 + , 9) will suffer 
of this SU (3) / breaking effect. In other words, |0 + , 9) can only be thought of as the leading term of the eigenfunction 
of Hcm-, without considering the correction from the next leading term caused by the strange quark content in + - 
tetraquarks. In a({ud}{ud}), there is no strange quark, so no such kind of corrections, hence \a) — |0 + , 9} a is suitable. 
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Figure 7: LHS of four interpolating currents of a+ and /o as functions of Borel mass squared with so=infinity in units of GeV 



In K,({ud}{ds}), there is one strange quark, its correction is relatively small, and the wavefunction |0 + ,9) K may be 
still valid to some extent. This is supported by numerical results. However, for fo{^{{ us }{us} + {ds}{ds})) or 

a+({us}{ds}), there are two strange quarks, the SU(3)f breaking effects is doubled. To these cases, one cannot insist 
the Jaffe's wavefunctions |/o) = |0 + ,9)/ o and \a + ) = |0 + ,9) a+ be still good enough to describe the non-perturbative 
QCD physics. Above all, we speculate that a legitimate SR analysis for /o and a+ should be based on the tetraquark's 
color-magnetic wavefunctions which are more precise, encoding the SU (3) ^-symmetry breaking effects. 




Figure 8: Mass of a+ and fo as function of Borel mass Mb with so being infinity. 
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IV. THE DIRECT INSTANTON CONTRIBUTION TO SUM RULE 



A. Analytic results 



In addition to the contribution of power type from the OPE expansion to the QCD SR, there are exponential 
contributions coming from direct instanton contributions. The direct instantion contributions originate from 't Hooft's 
instanton induced interaction [4J]. If the physics considered is relevant to two flavors, instanton effects induce a four- 
fermion interaction, as illustrated in Fig. [5] (usually called two-body single instanton contribution defined in [12]). 
In the framework of sum rule, this kind of instanton effect can be encoded in the quark propagator. Now the quark 
propagator has two terms, 

Sb = Sif> + S$~\ (71) 

S^fr corresponds to standard quark propagator (Eqs. ([53]) and ([54]) ) in Euclidean space, 5 , ^[ mst ' ) is related to instanton 
contribution and can be calculated by using the following formula in Euclidean space and regular gauge, 

S„f nSt) = M*> V)W»Q + ^)(Ur+r-U% b , (72) 

where 

r 2 

A q [x,y) = -i j^~2 ~ nx - z )(f)(y - zq) (73) 



167T 2 TO* 



and 



(j)(x-Z ) = - i , Ola/, - (74) 

l(x — Zq) + t \ ' 

Here r stands for the instanton size, zq for the center of the instanton. U represents the color orientation matrix 
of the instanton in SU(3) C and t+'~ are SU(2) C matrices. The effective mass of quark on the instanton vacuum is 
to* = m q — 2ir 2 r 2 (qq) /3 with current quark mass m q , here q € {u, d, s}. At the final stage, we multiply the result by 
a factor of two to take into account the anti-instanton effect and integrate over the color orientation and instanton 
size. When integrating over the instanton size, Shuryak's instanton liquid model [3l[ for QCD vacuum with density 
n r = n e ffd(r — r c ) has been used. 




Figure 9: The leading direct instanton contribution to the correlator, where "I" represents the instanton. 



With the definition Q 2 = —q 2 , the direct instanton contributions to the scalar nonet are listed below, corresponding 
to above two diagrams. Here, we only exhibit the contributions to cr-correlator, and the reader can find the results 
of other tetraquarks in appendix. We denote the total contributions from intanton and anti-instanton by "hist" . 
Recalling that the direct instanton contribution is possible only for different quark flavors, so in case of it, there is no 
direct three-body instanton contribution (from instanton induced six-fermion interaction). But to k, a+, /o, three- 
body instanton contribution might be important. However, in this paper, we only present the two-body instanton 
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contributions for these mesons, to capture the main physics 

T <T(inst) _ 156n eff r^(qq x2 



' - Trr fom (75) 



<r ] - ^4hiQ) + yf /owi. (76) 

7r s m*^ 187r 4 m* z 

<r , = _^£gl /6(Q) + I9y^y a ( } (?7) 

7r°m* z 187r 4 rriq 

n^ nst) = n^ nst) = 2 " e/ f% ( f )2 / (Q), (78) 

n# nst) = n^ lst) = 2nef J^ q)2 MQ), (79) 

H ^ " "7%^ /6(Q)+ /o(Q) ' (80) 

n .(inst) = n .(inst) = _^gi /6( Q), (81) 

II^ = n^ 1St) = - 2 -^lAf 6 (Q). (82) 



In above expressions, 



HQ) = J ^zo J dS 



iq-x 



fo(Q) = d z d x (83) 



The Borel transformation of Jq(Q) and fo(Q) are: 



o Jo 



B[MQ)} = — -r£- / d< / d y ^ r? — (X 2 + 5X d + 10X 



y 2 {i-y? 



-10X 5 + 5X 6 + X 7 ), 



B[fo(Q)\ = ^e- M «^/ 2 (Xo(Mir c 2 /2) + ^ 1 (Af|r c 2 /2)), (84) 
where we adopt the notations in paper [TJj], X = (1 — t)/t and K n {x) is the McDonald function. 

B. Numeric analysis of QCD sum rule with instanton effects 

To evaluate the direct instanton effects quantitatively, we make use of the following relation between the parameters 
of Shuryak instanton model [31j . 

Ueff 3 qe{u, d}, (85) 



m* 2 47r 2 r^ 



with 

r c = 1.6GeV _1 . (86) 
Considering the single instanton effects, the left hand sum rule becomes: 

n LHS (Q 2 ) = n OPE (Q 2 ) + n inst (Q 2 ). (87) 

After Borel transforming the both side of the QCD sum rule, we obtain the following relation 

B Mi U OPE (Q 2 ) + # M 2 n inst (Q 2 ) = 2nf x m* x e~ m x/ M * . (88) 
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In above expressions, 



B M 2_n 



OPE 



(Q 2 ) 



So 



e- s ' M Bp°™{s)ds, 



(89) 



where we have chosen a finite threshold to suppress the contribution from continuum. Utilizing the results in previous 
sections, the left hand sum rule can be performed for each possible interpolating current in (|43[) belonging to a certain 
meson. Then we can make use of the best current to fit the right hand sum rule to obtain the mass and residue. 
This approach was first suggested by (T^ |. In the following, for the sake of simplicity, we will only present a detailed 
analysis for a meson. For other mesons, the results are also exhibited. 

In Fig. [TQl we show the Borel transformed correlators IL{M%), including the instanton effects, at threshold value 
so=0.6 GeV 2 . From the Figure, we see that the instanton contributions are not always positive. To current Jf , they 
provide little negative contributions, and spoil the positivity of LHS obviously, when Borcl mass is small; to current 
J3 and J4 , instanton effects make the LHS rather negative, and this may be the usually called dangerous instanton 
contribution to sum rule [l2| ; only to current , the instanton effects improve the OPE calculation completely. This 
feature can be seen more clearly, if we notice that in Eqs. l|75p-(79): 



IT 



cr(inst) 
TT 



n 



cr(inst) 
PP 



n 



cr(inst) 
TP 



156n 



eff 



37r 4 mf 



HQ), 



32n 



7r 8 m* 2 



n 



r(inst) 
PT 



-HQ) 



19n e// r 4 (gq) 2 

187T 4 TO* 2 



2n 



eff 1 



\(qq)'< 



7r 4 m* 2 



HQ)- 



HQ), 



(90) 



In above expressions, the coefficients of fo(Q) are positive, while the coefficient of Jq(Q) is negative. After Borel 
transformation, HQ) and HQ) are J ust as m Eq. ([84]) . Numerically, B[f (Q)] is always positive, but B[f 6 (Q)] is 
always negative, so totally, the instanton contributions to the current J% are positive. From Fig. [lOl it is clear that 
the instanton contributions improve the convergence of current when Borel mass is small. At this moment, we 
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Figure 10: LHS of a including the instanton effects of four interpolating currents with so=0.6 GeV 2 



can use the numeric results associated with LHS of current J 1 ?, at threshold value 0.6 GeV 2 , to fit the RHS in single 
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Figure 11: The dashed and thick lines represent the left hand sum rule and right hand sum rule, respectively. To RHS, the 
mass and residue are presented in Table III. 

resonance approximation that is just the Eq. ([55)1 . as illustrated in Fig. [TTJ That choosing 0.6 GeV 2 as the value of 
threshold is inspired by previous OPE results. The fitted mass and residues are listed in Table 3: 

From the table, we notice that after adding the instanton contribution, the mass of a meson is still close to OPE 
result Eq. (|67l) . Then the instanton contribution is compatible with OPE results. It suggests that the physical mass 
of a depends weakly on the choice of QCD vacuum. 

Table III: Fitted masses and residues in single resonance approximation 



so (GeV 2 ) 


Ma (GeV) 


/ CT (HT 2 GeV) 


0.6 


0.72 


0.94 


1 


0.73 


0.93 



In the case of so = 0.6GeV 2 , considering a (±10%) variation of instanton size r c = 1.6±0.2, we find a corresponding 
variation of m a = 720^ 1 gQMeV and f a = 0.94l ° 7(10 _2 GeV). It seems like that the change of physical quantity lies 
within an acceptable range and the residue is more sensitive to the variation of instanton size compared with the mass. 
In [l5j], the authors discussed the meaning of the residue. In their notations, residue is defined as A 2 = 2ivf^Mf c . So 
we obtain a residue A 2 = 4 x 10~ 5 GeV 10 , which is larger than A 2 = 2 x 10~ 6 GeV 10 presented in According to the 
explanation of [15j |. large residue signifies the interpolating current operators have enough overlaps to the resonance 
states and the sum rule constructed with approximate OPE may contain enough information for the resonance to be 
extracted. So in our case, evaluating OPE up to dimension eight condensates seems reasonable. 

Finally, in order to investigate further the widths of the a meson states, it is necessary to find out three point 
correlation functions for er — > 7T7t, which has got out of the scope of this paper. 

As for other mesons, the current Ji still shows the best performance. The fitted masses and residues for k, a + and 
/o are presented in Table IV, V and VI in appendix , respectively. 

V. CONCLUSION AND DISCUSSION 

In this paper, we study the + nonet mesons as tetraquark states with interpolating currents induced from the 
color-magnetic wavefunction. This wavefunction is the eigenfunction of the effective color-magnetic Hamiltonian with 
the lowest eigenvalue, meaning that the state with this wavefunction is the most stable one and is most probable to 
be observed in experiments. Our approach can be recognized as constructing interpolating currents dynamically. We 
find that based on a current which is a proper mixture of the tensor and pseudoscalar contents, a good sum rule can 
be obtained. Our result can be perceived as a direct support to multiquark scenario described by the color-magnetic 
interaction, by means of QCD sum rule. 

In the SR calculations performed in this paper, we have taken into account the contributions from operators up to 
dimension d = 8 in the OPE. The results of SR analysis without instanton effects for + meson nonet {a, k, /o, a+} 
are : 
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1. a: In the SR analysis , a good Borel stability turns out in the region > 0.2 GeV 2 . Taking M\ f» 0.2 GeV 2 
and the threshold s ~ 0.6 GeV 2 , the largest pole contribution is 52% implying that a good SR analysis is 
achieved. Where we extract the mass of a (600 ± 75) MeV. 

2. k: A good sum rule was found when s = 0.9GeV 2 , M§ > 0.25GeV 2 . We obtain re mass (820 ± 80)MeV with 
pole contribution approaching 56%. 

3. fo and a + : to obtain a mass about 1 GeV by choosing the threshold and Borel mass, the pole contributions in 
SR are always around 10%. This indicates that the SR fails to analyze a+ and /o by using the interpolating 
currents deduced from the wavefunctions. We guess the reason is that in f (^({us}{us} + {ds}{ds})) or 

a+({us}{ds}), there are two strange quarks, so SU(3)f breaking effects are too strong to be negligible. This 
causes the Jaffe's wavefunctions |/o) = |0 + ,9)/ and \a + ) = |0 + ,9) a+ to miss some aspects of the fo- and a + - 
physics. We speculate that a legitimate SR analysis for them should be based on the tetraquark color-magnetic 
wavefunctions including the S'C/(3)/-breaking effects due to m c s > rh c . 

Proceed stepwise, we consider the direct instanton contribution. To the current Ji, the instanton effects are 
completely positive. Numerically, this positive effects improve the small Borel mass behavior of the Borel transformed 
correlator of current Ji. Meanwhile, adding instanton effects, the LHS gives a result compatible with OPE results. 

Finally, we go one step further and believe that the idea demonstrated in this paper also applies to 0~-g 3 <7 3 system. 
In [4l| . the authors have successfully extended Jaffe's method from q 2 q 2 to q 3 q 3 six-quark system (i.e., baryonium). 
One of the non-trivial results in [41] for baryonium is the existance of a counterpart of a. We denote this state by 
|0~,!y). Corresponding to Eq. ((HJ) for tetraquark, [4l[ shows 

H eff \0- ,l f ) = -82.533C|0- ,l f ). (91) 
In baryonium contents, its color-spin-flavor wavefunction can be expressed as: 

|0-,1/) = |1, 1/ ® = 0.591|(56 cs , 10 c ,4; 1/), (56 CS ,T0 C , 4; l f ), l c , 1; 1/ ® 1/) 

+0.807|(56 cs , 8 C , 2; 1 f ), (56 cs , 8 C , 2; 1 /), l c , 1; 1/ ® 1/), (92) 

where the notations in [4l| have been used. Like |er), |0~, If) has the largest mass defect among all the baryoniums. 
This implies that |0 _ , If), the lightest baryonium meson, may represent a stable physical state. Like Eq. ([5]), the 
mass of |0™, l_f) can be estimated roughly in the naive constituent quark model as follows 

m \o-.i f ) ~ \j)c^ l -~^ J ))\Q-,i f ) 

i i j 

(4 x 360MeV + 2 x 540MeV) - 82.533 x 4 X 20McV + 2 X 15MeV ^ 



« 1.007GeV. (93) 

We find that the mass of |0~, If) is close to that of 77' (960) [42^ . Furthermore, their quantum numbers are the same. 
So in the multiquark picture, we might identify |0~, If) as ry'(960) , or perceive 7/(960) as a baryonium or a Fermi- Yang 
meson [43[. Alternatively, there may be a large weight baryonium component in 77' (960). Usually, in the gg-picture, 
the mass of 7/ is attributed to U(1)a anomaly with non-trivial 9 vacuum in QCD [44J. However, that scenario has 
not excluded other schemes yet (e.g., see HH). In our case, a further examination to the conjecture on 77' in non- 
perturbative QCD should be meaningful. Since we have already known the color-magnetic wavefunction for |0 — , If), 
following the method presented in this paper, a SR analysis is straightforward. The result will be helpful to understand 
two interesting experimental measurements that may reveal the baryonium content of 77'. Those experiments are that: 

i) to measure the anomalous enhancement near the mass threshold in the pp invariant-mass spectrum from J/ip — > 
7PP reported by BES [46[ . 

ii) to observe resonance V(1835) in J/tp — > TK + Tr~rj' [47|]. In [46| the data fitting indicates that the enhancement 
is a S-wave Breit-Wigner resonance V(1835) [48|. It has been estimated that the decay branching fraction B(X — > 
pp) > 4% [49]. The decay mode of X — > pp is due to the tail effect of enhancement resonance of AT (1835) near the 
threshold of process J/ip — > jpp, therefore the fact of B(X — > pp) > 4% means the coupling between X and pp is 
very very strong. The most natural interpretation to this fact is that V(1835) is simply a bound state of p — p. 
Namely, V(1835) is a g 3 g 3 -baryonium molecular state [HO, IKH ■ In another hand, the major decay mode for A(1835) 
is A(1835) — * ir + n~r/' observed by BES [47| . It indicates that V(1835) is a molecular exciting state of meson 77' 
[5II ]. Consequently, the quark component of 77' should be same as A(1835), i.e.. 77' would be a 0~-baryonium meson, 
or a meson with large weight baryonium component. BES observations [46l . |47] provide evidence to this multiquark 
picture for 77' meson. 
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Appendix A 

1. Formulas of necessary spectral functions of k, a+ and fa 

For k ({ud}{sd}), since the current mass m s is much bigger than m u ,md, we can ignore terms proportional to 
m u , m<2 when listing the necessary spectral functions. Having done this, the length of formulas will be shortened, and 
the reader can have a clear impression about the structure of spectral functions. We will do the same thing for a + 
and /o- However, in numerical calculations, the contributions from the u,d quark mass terms have been taken into 
account. The spectral functions are the followings: 
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For a + ({us}{ds}) and f a (^{{su}{su} + {sd}{sd}), we only list the spectral functions for a + below. This is 
because in the widely adopted scheme Eq. ([53")) , u and d quark take the same value of current masses and condensates, 
which leads to a direct consequence that from the OPE calculation of the correlators of currents, we can not discern 
a+ and fo- In other words, to each kind interpolating current in Eq. (|43j) . the correlators of a+'s and the correlators 
of Jo's take the same expressions after completing the OPE calculation. 
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To convince the reader that our calculations are reliable, we make a comparison with the results of other authors. For 
example, 
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2. Instanton contribution to correlators of k, a + and /o 



We obtain the intanton contributions to k correlators as follows, 
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The instanton contributions to a + arc, 
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The instanton contributions to / are, 
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To check our results, we take the SU(3)f limit, which is m* = m* and (qq) = (ss). In this limit, the instanton 
contributions to n and a + are equal to each other. This is because that they all belong to the octet representation of 
SU(3)f. But this check is not suitable for / , since it comes from ideally mixing of the flavor singlet state with the 
isospin 1=0 component of flavor octet state . 



Table IV: Fitted masses and residues in single resonance approximation for k 



s (GeV 2 ) 


M K (GeV) 


/„( 10- 2 GeV) 


1 


0.72 


1.04 


1.5 


0.73 


1.00 



Table V: Fitted masses and residues in single resonance approximation for a+ 



so (GeV 2 ) 


M a+ (GeV) 


/ a+ (10- 2 GeV) 


1 


0.73 


0.89 


1.5 


0.73 


0.90 



Table VI: Fitted masses and residues in single resonance approximation for fo 



s (GeV 2 ) 


M /0 (GeV) 


/ /o (10- 2 GeV) 


1 


0.72 


0.90 


1.5 


0.72 


0.94 
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